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We study how the kinetic decoupling of dark matter within a minimal supersymmet¬ 
ric extension of the standard model, by adopting nine independent parameters (MSSM-9), 
could improve our knowledge of the properties of the dark matter protohalos. We show that 
the most probable neutralino mass regions, which satisfy the relic density and the Higgs 
mass contraints, are those with the lightest supersymmetric neutralino mass around 1 TeV 
and 3 TeV, corresponding to Higgsino-like and Wino-like neutralino, respectively. The ki¬ 
netic decoupling temperature in the MSSM-9 scenario leads to a most probable protohalo 
mass in a range of M p h ~ 10 _12 -1CU' M 0 . The part of the region closer to ~2 TeV gives 
also important contributions from the neutralino-stau co-annihilation, reducing the effective 
annihilation rate in the early Universe. We also study how the size of the smallest dark 
matter substructures correlates to experimental signatures, such as the spin-dependent and 
spin-independent scattering cross sections, relevant for direct detection of dark matter. Im¬ 
provements on the spin-independent sensitivity might reduce the most probable range of the 
protohalo mass between ~10~ 9 M 0 and ~10~' Mq, while the expected spin-dependent sen¬ 
sitivity provides weaker constraints. We show how the boost of the luminosity due to dark 
matter annihilation increases, depending on the protohalo mass. In the Higgsino case, the 
protohalo mass is lower than the canonical value often used in the literature (~10 _6 M Q ), 
while ( av ) does not deviate from (av) ~ 10 _26 cm 3 s _1 ; there is no significant enhancement 
of the luminosity. On the contrary, in the Wino case, the protohalo mass is even lighter, and 
(av) is two orders of magnitude larger; as its consequence, we see a substantial enhancement 
of the luminosity. 
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I. INTRODUCTION 

There is solid evidence that most matter in the Universe is in the form of non-baryonic dark 
matter (DM) [Th3]. From the theoretical point of view there are several particle physics theories 
which attempt to explain the yet unknown fundamental nature of DM. In the literature a plethora of 
DM candidates have been proposed (see, e.g., Ref. i). Depending on their masses and interaction 
cross sections with themselves or ordinary matter, they all exhibit a present day abundance in 
agreement with the DM density determined by Planck satellite, £l x h 2 = 0.1197T 0.0022 [6]. Among 
all particle physics candidates the most popular ones belong to the class of weakly interacting 
massive particles (WIMPs) [T, 7 ], since they are assumed to be stable and to have interactions 
with the standard model (SM) particles, giving a correct relic abundance as observed today. 

Although the SM describes the elementary particles and their interactions with great success, 
there are other good reasons, besides the need of a DM candidate, for expecting physics beyond 
the SM. One motivation is the so-called hierarchy problem. The mass of the Higgs boson acquires 
large quantum quadratic corrections proportional to the scale where the SM is valid. Assuming 
the SM is valid up to very high energy scales, the parameters in the theory have to be fine-tuned 
in order to keep the Higgs mass at an acceptable value of around 126 GeV. Since in the SM a 
symmetry that relates the various couplings does not exist, this situation is considered to be very 
unnatural (e.g., Refs. ISHEI). One of the best motivated scenarios introduced to solve this problem 
is supersymmetry (SUSY), with sparticle masses at the TeV scale. Although the first run of the 
Large Hadron Collider (LHC) placed important constraints to light sparticles, and a Higgs with 
126 GeV shifts the scale of SUSY to larger values requiring a certain amount of tuning (typically 
at 0(1%) for the MSSM, see e.g. [IT]). SUSY continues being a very attractive possibility. 

Another interesting feature of SUSY, mostly related to cosmology and the search for DM, is the 
existence of a conserved quantum number called R-parity, which assigns at each (super)partner 
of the SM particles R = — 1 while each ordinary particle is assigned R = +1. This quantum 
number implies that supersymmetric particles must be created or destroyed in pairs, and that the 
lightest supersymmetric particle (LSP) is absolutely stable, and hence DM candidates. In many 
supersymmetric extensions of SM, the lightest neutralino, a linear combination of the superpartners 
of the neutral gauge and Higgs bosons, is the favoured DM candidate. 

With the WIMP hypothesis, the abundance of DM originates from thermal decoupling in the 
early Universe. When the processes of pair-annihilation and pair-creation of WIMPs go out of 
chemical equilibrium due to the Hubble expansion, the resulting number density freezes out and 
remains constant per comoving volume until the present time. This chemical decoupling, however, 
does not signal the end of WIMP interactions with thermal plasma. There could still be elastic 
scattering processes with SM particles, which keep WIMPs in kinetic equilibrium until later time. 
When the rate for elastic scattering processes also falls below the Hubble expansion rate, WIMPs 
enter the epoch called kinetic decoupling. From this point on, WIMPs are decoupled from the 
thermal bath, and begin to free-stream. After this stage, first gravitationally bound DM structures 
begin to form, with the size set by the temperature of kinetic decoupling, related to a small-scale 
cutoff in the primordial power spectrum of density perturbations. Reference [12] calculated the the 
primordial power spectrum by including collisional damping and free-streaming of WIMPs, and 
showed that the free-streaming led to a cold DM (CDM) power spectrum with a cutoff around a 
scale corresponding to the Earth mass, 10 6 Mq (see also Refs. [THl - tTC] ). 

One of the most challenging goals today is to shed light on the nature of the small-scale cutoff in 
the primordial power spectrum of density perturbations, often dubbed with the name of protohalo 0 
Its properties are relevant for indirect DM searches. Indirect DM detection looks for signatures 


1 In the following, we use indistinctly “protohalos” or “subhalos” refering to protohalos, which are the smallest 
possible DM halos. 
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of DM annihilation, such as gamma-ray photons, from dense celestial environments, where the 
protohalo mass is a relevant quantity to determine the substructure “boost” factor. Direct detec¬ 
tion experiments of DM look for energy deposition in underground detectors caused by scattering 
interactions between target nuclei and WIMPs around us, giving valuable information about the 
scattering cross section, and through a correlation that we find in this study, they constrain the 
mass of the DM protohalos. 

Recently, Cornell and Profumo m studied scattering cross sections that are relevant for direct 
detection experiments and protohalo sizes in an MSSM context for the neutralino DM. They 
based their MSSM scan on 9 parameters defined at the electroweak scale. They found a strong 
correlation between the kinetic decoupling temperature and the spin-dependent (SD) cross section 
of neutralinos off nucleons. On the contrary, a weaker correlation was found in the case of the 
spin-independent (SI) neutralino-nucleon cross section. 

In the present paper, we do a forecast on the mass of the protohalos within a supersymmetric 
scenario by taking into account the latest data from all the relevant particle physics experiments 
as well as the relic density constraints. We perform our analyses within a Bayesian framework, 
by adopting 10 MSSM fundamental parameters defined at the gauge couplings unification scale, 
among which 9 of them we allow to vary after requiring the correct electroweak symmetry breaking. 

In the considered MSSM scenario, we find that the kinetic decoupling temperature leads to the 
protohalo mass most probably residing in a range of M p h ~ 10~ 12 -1CU 7 Mq. This large variation 
is due to the range of the kinetic decoupling temperature, T^d, since in the neutralino annihilation 
processes, both gauge bosons and fermions play a role, and these couplings reveal to be independent 
from one another. The range corresponds to two most probable posterior regions: Higgsino-like 
and Wino-like neutralinos, for which the most probable neutralino masses are around 1 TeV and 
3 TeV, respectively. 

In these most probable cases, we find that protohalo mass correlates with the both SD and SI 
scattering cross sections. We show that all Higgsino-like neutralino regions, where the probability is 
higher, such a scattering is dominantly spin-dependent. Therefore, any experimental measurement 
of the SD cross section will imply direct consequences on minimal protohalo rnassj^] 

We also show how future direct and indirect detection experiments can play an important 
role in constraining the (most probable) minimal protohalo mass down to rs_/ 10 9 M q and the 
expected value of the boost of the luminosity due to the annihilation of DM in those regions. 
Complementarity, we study how those predictions change in regions that are disfavored by the 
posterior probability density function (PDF) due to the large tuning, necessary to reproduce the 
experimental observables (including Mz)- 

This paper is organized as follows. We describe the supersymmetric model we adopt in Sec. [TTJ 
The role of the DM protohalo is discussed in Sec. m A brief explanation of the smallest DM proto¬ 
halo mass in Sec. Ill A The discussion of the most probable regions of the MSSM and interactions 
involved in the annihilation of neutralinos are presented in Sec. IIIB A profile likelihood map is 
discussed in Sec. Ill C We comment on the impact of the direct detection experiments on the mass 
of the protohalo in Sec. IV, and estimate of the boost of the luminosity due to the annihilation 
rate in a DM halo with substructures in Sec. [Vj We finally give our conclusions in Sec. VI 


II. THE MINIMAL SUPERSYMMETRIC STANDARD MODEL AFTER THE FIRST 

RUN OF THE LHC 

Despite the expectation around a potential discovery of light SUSY particles at the first run of 
the LHC, so far no signal of new physics has been found, which could be considered in tension with 

2 This is true if the scattering is mediated by a Z boson. The scattering could also be mediated by sleptons; in this 
case, we do not see such a correlation. 
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the ideas of natural SUSY. However, the relative large mass of the Higgs boson points to a heavier 
mass spectrum, suggesting that the lack of discovery of sparticles in the first run of the LHC is a 
consequence of the Higgs mass value. 

In the MSSM, a lightest Higgs boson of around 126 GeV implies a range of MgusY between 10 3 
GeV and ~ 3-10 4 5 GeVj^Jwhere ALsusy represents the scale at which SUSY particles decouple from 
the SM (for details, see |18U19| ). Hence, within the MSSM framework the Higgs mass is in tension 
with naturalness of the electroweak symmetry breaking, requiring a typical tuning of 0(1)%, see, 
e-g-, HU. This tension is relaxed going beyond the MSSM [20H24] . Moreover, since stops give the 
most important contribution to the Higgs mass, the allowed range of AfgusY could be written as a 
constraint to the stop sector, where typically stop masses should be larger than ~3 TeV, unless its 
mixing parameter reaches its maximal value [25] , leaving basically the rest of the SUSY spectrum 
unconstrained 0 

On the other hand, one of the beautiful aspects of SUSY is the apparent unification of gauge 
couplings in the MSSM, because it gives a strong hint in favor of grand unified theories suggesting, 
as well, that we know how the renormalization group equations (RGE) behave up to the gauge 
coupling unification scale, A/gut 13 Taking SUSY parameters at ALgut leads to implicit relations 
between sparticle masses, in particular the average of stop masses at the scale of 1 TeV for tan /3 = 

10, mf , written as a function of the soft parameters at ALgut reads, HU: 

£ 1,2 

m? 2 ~ ( 2.972 Aff + 0.339 m 2 ^ + 0.305 m 2 ^ + 0.091 Aff - 0.154 m 2 Hu - 0.052 A 2 (1) 
+ 0.017M? ... ) + m 2 t , 

where M\, M 2 and M 3 are the bino, wino and gluinos soft mass terms, respectively, and 
are the third generation of squark soft masses, and rri h. u is the H u soft mass. Equation 0 shows 
that large stop masses imply large gluino mass (Mg ~ 2.22M 3 ), unless the soft mass terms of the 
third generation squarks are very large, which leads to a scenario like split SUSY HU. 

Regarding naturalness, the largest tuning required to get the correct electroweak symmetry 
breaking is applied on the g parameter. From the minimization of the Higgs potential one obtains 

^ M | = ( 1.62 Aff - 0.64 rn 2 Hu + 0.37 m 2 Q3 + 0.29 rnfj 3 - 0.29 A t M 3 - 0.20 Aff (2) 

+0.14 M 2 M 3 + 0.11 A 2 t +... )-g 2 , 

where this expression is valid at a scale of 1 TeV for tan/3 = 10 HU- As in Eq. M 3 is the 
responsible for the larger contribution. The current gluino mass bound from ATLAS and CMS 
[251 29]. rrig > 1.33 TeV (assuming 100% decay to qq\\ and a mass difference between g and °f 
at least 200GeV), is that more stringent for naturalness. From Eqs. 0 and Q, we could also see 
that naturalness and Higgs constraints affect mainly the gluino and squarks sector. On the other 
hand, sleptons, Binos and Winos are basically unconstraint. 

In a more general framework, where the MSSM is parameterized at EW symmetry breaking 
scale, the pMSSM, the Higgs mass measurements constrain mainly the stop sector, leaving the rest 
of the spectrum effectively unconstrained. In this case, the main constraints for sparticle masses 
come from LHC limits and B-physics (see, e.g., Refs. [30H33] 1. 


3 This range is valid for relatively large values of tan /3. 

4 Notice that Ref. El re-examined the natural SUSY scenarios, and showed that light stop masses (closer to its 
lower limit after imposing the Higgs mass) are not really a generic requirement of natural SUSY scenarios. 

5 I 11 gravity, mediated SUSY breaking scenarios conditions are set at Mpi an k- A popular approximation is to start the 
RGE running from Mgut instead of Mpi an k. For some particular scenarios, this approximation is not necessarily 
correct [26] . 
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Besides the tuning associated to the EW symmetry breaking, there is also a tuning associated 
to the requirement of having a good DM candidate. Refs {331 SSJ study the fine tuning required 
to obtain the correct DM relic density. In particular, |35] shows that the region of 1 TeV, corre¬ 
sponding to the lightest Higgsino-like neutralino, requires very smallest tuning. Typically, regions 
where the correct annihilation cross section is dominated by resonances or sfermion-neutralino 
co-annihilations require a large tuning. 

To study the MSSM parameter space we perform a Bayesian analysis. One of the interesting 
aspects of this approach is that it is possible to take into account naturalness arguments [36]. A 
fine-tuning associated to the electroweak symmetry breaking is included when we takes the mass 
of the Z boson in the same foot as rest of the experimental data. Effectively, after requiring the 
correct electroweak symmetry breaking, the posterior PDF appears to be proportional to a term 
that penalizes regions with a large fine-tuning, independently of the choice of the prior probability. 
Interestingly, this term is inversely proportional to the Barbieri-Giudice fine-tuning parameter Ell- 
More specifically the EW fine tuning penalization appears as a Jacobian factor that arise from the 
change of variables {gi, yi, g, B} -» {a*, rrif, Mz, tan/3} evaluated at the measured value of Mz, 
where gi and y* are the gauge and Yukawa couplings respectively, B is the bilinear Higgs coupling, 
and /x is the Higgs mass term in the superpotential defined at the SUSY breaking scale. This 
Jacobian factor is completely independent of the choice of parameters and is not based in a specific 
definition of fine-tuning. In the same way, a fine-tuning penalization associated to all the other 
experimental observables is included. Motivated by the fact that this definition does not involve 
prejudices, Ref. |38| came up with the idea of using this covariant matrix to penalize regions with 
large fine-tuning in a x 2 analysis. 

In our analysis we assume gravity mediated SUSY breaking and parameterize the MSSM with 
10 fundamental parameters defined at the unification scale of the gauge couplings as well as SM 
parameters. We also assume unification and universality conditions for the squark masses, slepton 
masses and trilinear terms. The set of 10 parameters is: 

{gi,Vi, Mi, M 2 , M 3 ,rn 0 q ,mo l ,m H , A s 0 q , A s 0 l ,g,B) , (3) 

where we added, as well, the gauge and Yukawa couplings, y* and y u respectively. Mi, M 2 , M 3 
are the gaugino masses, m s 0 q , rriQ and mu are the soft squark, slepton and Higgs masses, A s 0 q and 
K are the squarks and slepton trilinear couplings, B is the bilinear Higgs coupling, and g is the 
Higgs mass term in the superpotential. 

Using a more convenient parameterization, the effective set of parameters reads: 

{s, Mi, M 2 , M 3 , m s 0 q , mo, m H , A s 0 q , A s 0 l , tan/3, sign(y)}, (4) 

where s stands for SM parameters described in Table [T] and, without loss of generality, the sign of 
/j is fixed to +1, allowing Mj to have positive and negative values. In such a way we cover regions 
with relative phases between g and Mj. 

Let us comment about how strong the predictions of the scenario we consider are with respect 
to the most general MSSM. In our approach we assume that SUSY was broken at gauge-coupling 
unification scale. Although this assumption is reasonable in gravity-mediated SUSY breaking 
scenarios, it is not the only possibility; for example, in gauge mediated scenarios it can happen 
in principle at any scale. Moreover, the consequence of this assumption depends on the freedom 
we give to the soft parameters. Imposing universality condition (squark and slepton squared-mass 
matrices proportional to the 3x3 identity matrix) and unification condition (right sfermion masses 
equal to left sfermion masses and = m# d ) implies a specific mass hierarchy for squarks and 
sleptons which could be ameliorated if the SUSY breaking scale is smaller. Hence, the regions 
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Gaussian prior Range scanned 

ref. 

M t [GeV] 

173.2 ±0.9 (167.0,178.2) 

m 

m b (m b ) MS 

[GeV] 4.20 ± 0.07 (3.92, 4.48) 

[43] 

[aem(Mz) MS ] _1 127.955 ± 0.030 (127.835, 128.075) [43] 

a s {M z ) MS 

0.1176 ± 0.0020 (0.1096, 0.1256) 

m 


TABLE I. Nuisance parameters adopted in the scan. 

of parameters we are missing in using this parameterization of the MSSM are the ones with any 
possible hierarchies of sfermions masses. In our case, t\ is always the lightest stop and f\ the lightest 
slepton. On the other hand, the universality condition is supported by the strong constraints from 
FCNC process. 

Using a more general parameterization at EW symmetry breaking scale, the pMSSM, the spar- 
ticle masses do not feel the impact of the renormalization group equations J^] the correlation between 
the parameters disappear and the choice of the prior will most likely dominates the results not 
allowing us to make conclusions about the most probable region. On the other hand, Bayesian 
analysis has been performed in the pMSSM from different perspectives, to be able to identify 
which are the parameters that are directly constrained by the experimental information, that can 
be checked by looking at the prior dependency in each parameter (see, e.g., Refs. [39HH]). 

To perform the analysis, we follow the lines described in Ref. [55], where two different priors are 
considered: standard log priors (S-log prior), which takes a log prior for each parameter indepen¬ 
dently, and improved log priors (I-log prior), which assumes a common origin for the soft-masses, 
as expected from SUSY breaking mechanisms. The range of the parameters in our scan varies from 
10 GeV to 10 6 * 8 GeV. Although both of the considered priors are based on logarithmic space, they 
are quite different from one another; S-log prior, for example, favors large splittings between the 
parameters, while I-log priors assume a common origin for the soft parameters. For a more detailed 
discussion about the priors see Sec. 3.3 of Ref. [45]. Notice that following this approach, which 
takes naturalness arguments into account, we are able to explore a large range of the parameters 
and get a consistent result. Previous Bayesian analyses followed a different approach finding prior 
dependency in their results, showing that not including M z as a experimental observable, and 
therefore not taking into account EW fine-tuning, it is not possible to conclude about the most 
probable region, for example in the CMSSM. 

The experimental data considered in our analysis is described in Table [IT] where we include 
electroweak precision measurements [46], B-physics observables the Higgs mass [561 Ef], 

and constraints on the WIMP-nucleon scattering cross-section by XENON-100 |58| . In addition, 
we include the measured relic density according to Planck results [59] because we assume a scenario 
with a single DM component which is produced thermally in the early Universe^ 

For the numerical analysis we use SuperBayeS-v2.0, a publicly available package that include 
MultiNest m [68] nested sampling algorithm, Softsusy [65] for the computation of the mass 
spectrum, micrOmegas [70] for the computation of the relic density, DarkSusy [7T] for the compu¬ 
tation of direcl|^]and indirect detection observable, SusyBSG [73] and Superiso [74] for B-physics 
observable. 

6 However, the universality condition is somehow taken into account in the pMSSM, when setting first and second 
generation sfermion masses equal. 

' The updated values for B decays are, for example, BR(B — > sy) = (3.43 ± 0.22 ± 0.07) x 10 -4 [52] (see also 
Refs. [531 l54j ) and BR(B — > a + n~) = 2.8127 x 10~ 9 [55]. All these measurements are still in agreement (within 
uncertainties) with the values that we adopted in our analysis, and therefore, their impact would not be large. 

8 In our analysis, we assume that 100% of dark matter consists of the neutralino. If there is other dark matter 
components, we need to regard the measurement of the dark matter density determined by Planck satellite as an 
upper limit, and follow some scaling ansaz studied in, e.g., |45ll60l . This is however beyond the scope of this paper. 

9 For the contribution of the light quarks to the nucleon form factors, concerning the spin-independent WIMP- 
nucleon cross section, we have adopted the values /t„ = 0.02698, frd = 0.03906 and /t s = 0.36 derived 
experimentally from measurements of the pion-nucleon sigma term. 
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Observable 

Mean value 

Uncertainties 
a (exper.) r (theor.) 

Ref. 

M w [GeV] 

80.399 

0.023 

0.015 

m 

sin 2 9 eff 

0.23153 

0.00016 

0.00015 

ED 

BR (B X s -f) x 10 4 

3.55 

0.26 

0.30 

m 

Mb s 

BR(B u ^tis) 

1.04 

0.11 

- 

HU 

1.63 

0.54 

- 

m 

A 0 _ x 10 2 

3.1 

2.3 

- 

[M] 

BR {B^Dtv) „ lf)2 

BR (B^Dev) A lu 

41.6 

12.8 

3.5 

m 

R/23 

0.999 

0.007 

- 

p] 

BR(D S —» tv) x 10 2 

5.38 

0.32 

0.2 

m 

BR(D S —> fj,v) x 10 3 

5.81 

0.43 

0.2 

m 

BR(L> jiv) x 10 4 

3.82 

0.33 

0.2 

m 

n x h 2 

0.1196 

0.0031 

0.012 

m 

m h [GeV] 

125.66 

0.41 

2.0 

m 

BR(R S -> /i + p“) 

3.2 x 10" 9 

1.5 x 10" 9 

10% 

ED 


Limit (95% CL ) 

r (theor.) 

Ref. 

Sparticle masses 
m x - af N 

As in Table 4 of Ref. [65]. 
XENON100 2012 limits (224.6 x 34 kg days) 

|MJ 


TABLE II. Observables used for the computation of the likelihood function. For each quantity we use 
a likelihood function with mean /i and standard deviation s = V a 2 + r 2 , where a is the experimental 
uncertainty and r represents our estimate of the theoretical uncertainty. Lower part: observables for which, 
at the moment, only limits exist. The explicit form of the likelihood function is given in ref. [65]. In 
particular, in order to include an appropriate theoretical uncertainty in the observables, the likelihood 
contains a smearing out of experimental errors and limits. 

For the Wino-like and Higgsino-like LSP cases, the Sommerfeld enhancement^ of the primordial 
and present day neutralino annihilation has been included, following the lines of Refs. C3M, 
using DarkSE EBEa, which is a package for DarkSusy. We created a grid in the plane and 

performed interpolations to correct the values of the relic density and the present day neutralino 
annihilation within SuperBayeS interface. 


III. DARK MATTER PROTOHALOS IN THE MSSM 


A. The smallest mass of the protohalo 


WIMP interactions with the plasma in the early Universe produce damping of the power spec¬ 
trum before and after the kinetic decoupling. Before kinetic decoupling, WIMPs behave as fluid 
tightly coupled to the plasma. Interactions produce shear viscosity in the WIMP fluid causing the 
density perturbations in the WIMP fluid to oscillate acoustically in the heat bath mm- The 
damping scale set by acoustic oscillations is given by the DM mass enclosed in the horizon at this 
epoch, i.e. the size of the horizon at kinetic decoupling [55j : 


M n 


47r Px( T kd) 

3 HHT kd ) 


= 3.4 x 10“ 6 M 0 


I T kd gj ff 

l 50MeV 


(5) 


10 The Sommerfeld enhancement [75] is a nonrelativistic effect that depends on three quantities: the neutralino mass, 
the difference in mass between the neutralino and the next to the lightest particle, and the size of the coupling 
among them. In this context, Sommerfeld enhancement of the annihilation cross sections can significantly shift 
the neutralino mass consistent with the experimental Qcdmh 2 value l76j . 
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where g e ff is the number of effective degrees of freedom in the early Universe and p x is the DM 
density, both evaluated at the temperature of kinetic decoupling, T^. 

After kinetic decoupling, WIMP interactions give a free-streaming scale which induces a damp¬ 
ing of density perturbations below a scale characterized by a (comoving) free-streaming wavenum¬ 
ber, kf s [T2'l [M]. Therefore, if we have perturbations contained in a sphere of radius n/kf s , we 
have the minimal mass a DM protohalo, and then the mass of the smallest protohalo allowed by 
free-streaming is [83]: 


1 + ln(5^T M /50MeV)/19.1 
(m x /100 GeV) 1 /2 ffe 1 /J(r M /50 MeV)V2 

The mechanisms of collisional damping and free-streaming of WIMPs lead to a cutoff in the CDM 
power spectrum, from which the typical scale for the first haloes in the hierachical picture of 
structure formation is set. The canonical value for the mass of the DM protohalos is related to the 
nature of the DM particle. The SUSY prediction for the size of the DM protohalos falls in a range 
from 10~ 11 to 10 -3 Mq [83]. It is not clear if these first and smallest halos survive until today, 
since they can be destroyed either in the process of merging or by star formation. According to 
Refs. [[MI f85 j|. the first halos lose their mass during structure formation, but survive until today 
with their inner density still intact. 

For a typical WIMP one finds that the chemical-decoupling temperature is given by T c d ~ m/25, 
where m is the WIMP mass, and the annihilation cross section by (av) ~ 10 _26 cm 3 s _1 in order 
to obtain the CDM relic density as observed today. On the other hand, the kinetic-decoupling 
temperature and, therefore, the minimal protohalo mass are not well constrained for WIMPs. 
A reference value of the minimal protohalo mass for SUSY candidates is 10 6 Mq, which was 
computed by assuming a Bino-neutralino scattering with the SM particles through sfermions with 
a mass of around twice the neutralino mass HUES [86] [87]. The chosen nature of the neutralino 
and the particular relation between the sfermions and Bino neutralino was well motivated by 
constrained SUSY extention of the SM (CMSSM), where the typically light neutralinos (lighter 
than 1 TeV) that are able to reproduce the correct relic density are mostly Bino that, efficiently, 
annihilate through sfermions in the early Universe. 

Even though M p h ~ 10 -6 Mq was a good estimate of the value of the smallest mass of the 
DM protohalos for a “well motivated” neutralino, it is not a strong prediction for a general neu¬ 
tralino DM. As described in [88], there are several ways to get a well tempered neutralino. The 
Bino-neutralino, that annihilates through sfermions, is one of those possibilities. Reference [83] 
performed a general study of the T c d and T^d for the MSSM neutralino, as expected m/T c d ~ 25, 
while m/Tkd has a range of variation of almost four orders of magnitude, leading to a range in 
Mph ~ 1O- 12 -1O- 3 M 0 . The reason of the big range of T^d is that the interactions involved in 
the annihilation of neutralinos, that are constrained by the relic density, are not necessarily the 
relevant for the last scattering of neutralinos with the plasma, and therefore the relic density does 
not constrain it. For example, in the case of Wino-like or Higgsino-like neutralinos, the annihilation 
products are mainly gauge bosons, whose interactions involve different couplings with respect to 
the ones involved in the neutralino-fermion scattering. 

It is important to mention that the computation of the kinetic decoupling temperature, and 
hence, the smallest protohalo mass becomes more complicated when the decoupling occurs close to 
the Quantum Chromodynamics (QCD) phase transition. As mentioned above, the computation of 
T^d and M p h was performed with DarkSusy, following the lines described in Ref. [83]. For the case of 
two light (u, d ) and one massive s quarks, the critical temperature is assumed to be T c = 154 MeV. 
The plasma is described including three quarks and gluons for a temperature T > 4T C . Therefore, 


47r 


7r 




M fs ~ ~7rPx ( t~ ) = 2.9 X 10 -6 M e 


o 
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in the following analysis, for the regions where the kinetic decoupling temperature lies between T c 
and 4 T c (154 MeV < < 616 MeV), T^d will represent an upper bound while M p h a lower bound. 


B. The most probable regions 

The determination of the smallest mass of the DM protohalo for the most probable regions of 
the MSSM is of great interest for the study of both direct and indirect detection of DM. 




FIG. 1. The two dimensonal joint posterior probability density function for the temperature of kinetic 
decoupling, Tkd, versus the neutralino mass (left panel), and for the protohalo mass, M p h, versus the 
neutralino mass (right panel). The region with higher probability density corresponds to a Higgsino DM 
candidate, while in the second region the DM candidate is a Wino. 

In Fig. [l] we show the two dimensional joint posterior PDF for the temperature of kinetic 
decoupling, T^d, and for protohalo mass, M p h, against the neutralino mass. The contours represent 
intervals at 68% and 95% credible regions. The two most probable regions are around ~1 TeV and 
~3 TeV and correspond to Higgsino-like and Wino-like neutralino, respectively. 

As discussed in pR>j, both Higgs mass measurement and relic density constraint are the main 
responsible for the shift of the preferred regions towards higher masses. We would like to stress that 
the credibility intervals represent the most probable region assuming the model that we consider 
is correct. In other words, the credibility regions show the relative probability density within the 
model. Points outside the contours are disfavored because they have worse likelihood and/or they 
require a large tuning to reproduce the experimental data. Concerning the prior dependence of 
our Bayesian analysis, we checked the stability of our results by using two different priors (I-log 
and S-log priors), finding that the result is basically the same; it means that our result is prior 
independent. 

Let us describe in more details the two most probable regions. The region around 1 TeV 
corresponds to a Higgsino-like neutralino, whose annihilation cross section is driven by its Higgsino 
component (the main annihilation processes are those of a pure Higgsino-neutralino). On the 
other hand, for the scattering cross section, the small component of Wino and Bino plays a crucial 
role. The reason is the following. Assuming that sfermions are decoupled, the tree level SD 
scattering of Higgsino-like neutralino with fermions is mediated by the Z boson. In the limit 
of pure Higgsino-neutralino, H u and Hd are degenerate, and since they have opposite quantum 
numbers, their contributions cancel. But then, when the gaugino masses are not decoupled, the 







10 


H u and H^ composition of the lightest neutralino is not the same, and the cancellation does not 
occur. Regarding the tree level SI scattering]^] the interaction is mediated by the Higgs boson, and 
as it interacts with the neutralinos via a Higgsino-Bino(Wino)-Higgs coupling, a nonzero gaugino 
component is necessary in order to have a nonzero tree level contribution. In this region sfermions 
are not necessary heavy enough to be considered decoupled. However, since Higgsino-sfermion- 
fermion interaction is proportional to the Yukawa coupling]^] these contributions are negligible. 

The region around 3 TeV corresponds to Wino-like neutralino, where the most important an¬ 
nihilation interactions are those of the pure Wino neutralino. The part of the region closer to ~ 
2.5 TeV has also important contributions from the neutralino-stau co-annihilationp] reducing the 
effective annihilation rate of neutralinos in the early Universe and, therefore, decreasing the value 
of the neutralino mass to obtain the correct relic density, that for the case of pure Wino is ~3 TeV. 

As in the Higgsino-like neutralino case, the tree level SD neutralino-sfermion scattering cross 
section receives an important contribution from the Z boson, which is the mediator of this inter¬ 
action; while the tree level SI neutralino-fermion scattering cross section from a Higgs. In both 
cases, a non-negligible component of Bino or Higgsino is needed to have a tree-level contribution to 
these processes, since W°~W°-Z and W°-W°-h interactions do not exist. In addition, sfermions 
give an important contribution to the neutralino-fermions scattering cross sections, in particular 
for Wino-neutralinos with mass ~ 2.5 TeV. As we commented above, in this region staus are close 
in mass to the lightest neutralino, and selectrons and smuons are light enough to give a sizeable 
contribution to the scattering cross section. 

C. Profile likelihood maps 

In the previous section we showed that the most probable neutralino mass regions are those 
around 1 TeV and 3 TeV. We once again underline that this result is based on the relative proba¬ 
bility density between the regions of the model. It does not imply that there are no valid points in 
the region of lighter neutralinos, i.e. in the intermediate region between 1 TeV and 3 TeV. In order 
to have points with good likelihood outside the 95% credibility region showed in Fig. [lj a (larger) 
fine-tuning which reproduces both the experimental data and the correct electroweak symmetry 
breaking is required. 

In this subsection we study models that reproduce all the observables within 2cr confidence 
level. To this end, we performed a new exploration by requiring a non-negligible Bino component 
for the lightest neutralino. In this way we completed our previous exploration related on the study 
of the Higgsino-like and Wino-like neutralinos, including all the different neutralino natures. We 
included some of the latest ATLAS bounds on sparticle masses based on simplified models detailed 
in Table [TTTj To apply the simplified model limits, we use the production cross sections published by 
LHC SUSY Cross Section Working Group [90j, which performs an interpolation routine for gluino, 
squark and neutralino-chargino production. Slepton production cross section has been computed 
using PYTHIA 8 |91L [92l . We also include the overall signal strength of the Higgs measured by 
ATLAS [93]. For the computation of the branching ratios we used SUSY-HIT [91] . 

Figure [2] shows points that reproduce the experimental data within 2 a confidence level. We 
show the lightest neutralino mass as a function of the kinetic decoupling temperature, T^i (top), 
and the protohalo mass, Af p h (bottom). Let us describe the mass spectrum. The characteristics 

11 We still assume that sfermions are decoupled. 

12 We remind that at temperatures of the order of MeV, when the kinetic decoupling occurs, the population of third 
generation of fermions is very small. 

13 Sometimes solving the Boltzmann equation for the evolution of the neutralino number density to obtain the correct 
relic abundance of DM requires additional considerations; degeneracies in mass between the lightest neutralino 
and the next to the lightest one, or the presence of thresholds and resonances in the annihilation cross section 
may be relevant (see, e.g., the review [55]). In particular, when the lightest neutralino is close in mass to a heavier 
neutralino, the relic abundance is determined both by its annihilation cross section and by co-annihilation with 
this heavier partner that, then, decays into the lightest one. Co-annihilations may also occur with squarks, when 
they happen to be very close in mass to the lightest neutralino. 
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Topology 

Luminosity 

Reference 

Production 

Decay 

Comment 

titi 

h bW^xl 

m; <C m + 

20.3 fb' 1 

m 

tih 

ti -> tx i 

all hadronic 

20.1 fb _1 

m 

bib t 

bi ->■ bx\ 


20.3 fb _1 

m 

99 

9 -> bb\ i 

rriq 3> m,g 

20.1 fb _1 

M 

99 

9 -> tt \i 

rriq Trig. 0 leptons + 3 b-jets channel 

20.1 fb" 1 

m 

99 

9 99X i 

rriq S> mg 

20.3 fb" 1 

m 

99 

9 bt X f 

mg mg, to, ± — to, o = 2 GeV 

20.1 fb _1 

m 

99 

9 9Xi 

TOg » mg 

20.3 fb _1 

m 

Xtx ° 2 

W^ X °iZ^X°i 

X 

II 

X 

too 

20.3 fb" 1 

Ei 


II 


20.3 fb _1 

m 


Ir i ± X? 


20.3 fb" 1 

m 

^LR^LR 

%R ->• i ± X? 


20.3 fb" 1 

m 


TABLE III. Simplified models exclusion limits we have included in our analysis. 


of the electroweakino sector are set mainly by the fact that an efficient neutralino annihilation is 
needed to reproduce the correct relic density. To identify regions where the lightest neutralino co¬ 
annihilates with sfermions in the early Universe, in both left panels we highlight points that satisfy 
a criterion based on the mass difference between the lightest neutralino and the lightest stau (green 
points), and between the lightest neutralino and the lightest stops (blue points). To select those 
points we have required a maximal relative mass difference, A (mj — m x o), of 5% and a maximal 
absolute mass difference of 5 GeV which are imposed for neutralino masses above and below 100 
GeV, respectively. The gray band of the top-right panel shows the range of temperatures where 
the QCD phase transition occurs, from the critical temperature to four times this one, where the 
value of Tm represents an upper bound. Those points with a T^ around the QCD phase transition 
are represented with lighter colors in the m x o-M p h plane in the bottom-left panel, where in this 
case the value of M p h represents a lower bound. 

Most of the points showed in Fig. [2] have a neutralino quasi-degenerate with another sparticle. 
Higgsino-like and Wino-like lightest neutralinos are quasi-degenerated with the lightest chargino, 
guaranteeing both a very efficient annihilation of neutralinos and co-annihilation with charginos, 
and selecting rather heavy neutralino masses. Neutralinos with a dominant Higgsino or Wino 
component cover the mass region of m x o > 1 TeVp*| As we commented in the previous section, 
for pure-Higgsino and pure-Wino neutralino the relic density constraint fixes the mass to ~1 TeV 
and ~3 TeV, respectively. As a result of our scan we have identified different mixed states lying 
between these regions: Higgsino-Wino neutralinos, Higgsino-like and Wino-like neutralino that 
co-annihilate with staus or stops, and Wino-Bino neutralinos and Higgsino-like neutralino with a 
mass equal to half of the mass of the pseudoscalar. 

Some of the points with m v o slightly below 1 TeV are Higgsino-Bino neutralino. This region is 
strongly constrained by direct detection experiments like XenonlOO and LUX. However, there are 
some blind spots for /j, and M\ with opposite relative sign, as explained in detail in Refs. (31] [35]. 

Points with m v o < 1 TeV have a lightest Bino-like neutralino. For 100 GeV < m. o < 600 GeV 

X-l _ . _ ^J_ 

it is possible to distinguish two groups of points in the bottom-left panel of Fig. [2] The first 
group has smaller T^, ranging from ~10 MeV to ~100 MeV and is basically aligned to the stau 
co-anhihilation region. For these points sleptons are light, and the correct neutralino abundance 
was reached by slepton-neutralino co-annihilation in the early Universe. Charginos and heavier 
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Assuming DM is made of several species, the relic density constraint becomes an upper bound, allowing to have 
lighter Higgsino-like and Wino-like neutralinos. 
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100 1000 100 1000 
m v o (GeV) m y o (GeV) 

X\ 

FIG. 2. Lightest neutralino mass versus Tkd (top panel) and M p h for points that reproduce all the experi¬ 
mental observables within 2er confidence level. 


neutralinos are typically much heavier. For the second group of points with larger Tkd varying 
from ~100 MeV to ~3 GeV we checked that the lightest (Bino-like) neutralino is quasi degenerated 
with both the lightest Wino-like chargino and the second lightest neutralino, guaranteeing the 
neutralino annihilation. Top-right panel of Fig. [2] shows that these two regions are not completely 
disconnected. For example, for A(m ; ~ — m x o) ~ 0.5 (meaning mj ~ 3m x o) sleptons also play a role 
in the annihilation processes. 

The region 600 GeV < m x o < 1 TeV has similar characteristics, but in this case the two 
regions, that one with light sfermions and the other one with light chargino, have a large overlap 
for 30 MeV < Tkd 1$ 500 MeV. 


Last but not least, we find that there are very few points for the Higgs and Z resonance regions. 
These two regions require a very large tuning, and therefore, they are very difficult to explore when 
requiring boundary conditions at GUT scale. 

To understand the dominant process of neutralino-SM scattering in the regions we described 
above, in the top-right panel of Fig.[2]we show the relative mass difference between the lightest first 
and second generation of sleptons and the lightest neutralino, A(m ; - — m x o), while in the bottom- 
right panel we show the gaugino fraction 15 These plots show, for all gaugino-like neutralinos 


lo The lightest neutralino is a linear combination of the superpartners of the gauge and Higgs field: \i = NuB + 
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(Bino-like or Wino-like), a clear correlation between the lightest neutralino mass and the kinetic 
decoupling temperature for a fixed value of A (mj — m x o). Higgsino-like neutralinos around 1 TeV 
do not show a correlation for specific sleptons masses, as we comment in the previous section; its 
interaction with sfermions is proportional to the Yukawa coupling, and it is therefore negligible for 
the first and the second generation of sleptons. In the Higgsino-like case the dominant interaction 
is the one mediated by the Z-boson, as in the case of Bino-like and Wino-like neutralinos when 
sfermions are decoupled. 

As we commented in section |TTJ we assume universality and unification of squarks and slepton 
masses. These conditions imply that the t\ and f\ are the lightest squark and slepton, respectively, 
which is the reason why we only find neutralino-stop and neutralino-stau co-annihilation regions 
in our analysis. In more general scenarios where sfermions masses do not unify, the possibility of 
having co-annihilation with any sfermion is open, since any of them could be the next-to-LSP. If 
the lightest neutralino is Bino-like and the first or second generation sfermions are close enough 
in mass to the lightest neutralino to guarantee a large enough effective annihilation in the early 
universe then the dominant interaction in the scattering between the lightest neutralino and the SM 
particles will be the same interaction (neutralino-fermion-sfermion), producing strong correlation 
between the mass of the lightest neutralino and M p h- 

Another consequence of universality and unification is that the first and second generation of 
squarks are in general very heavy (due to the Higgs mass constraint to the stop sector), having, in 
most of the cases, a negligible contribution to the neutralino annihilation and neutralino scattering 
with the SM particles in the early universe. Without this assumption the most important constraint 
to squark masses will come from LHC bounds and direct detection experiments, allowing smaller 
masses. Due to the strong lower bounds on first and second generation squarks masses coming from 
LHC [ 101) , one will expect that sleptons will still give the dominant contribution to the neutralino 
annihilation and neutralino scattering with the SM particles in the early universe for neutralinos 
lighter than 300 GeV. However, for neutralino masses larger than 300 GeV, contributions from 
first and second generation squarks could be sizeable. 

Interestingly, the cases that set the smaller value of M p h, when sleptons are very close in mass to 
the mass of the lightest neutralino, and larger value of M p h, when sleptons are decoupled and the 
scattering is mediated by Z-boson, are covered in our analysis. On the other hand, the consequence 
on LHC, direct detection and indirect detection could be different, as we will discuss in the next 
section. 

The understanding of the interactions that play a relevant role in the annihilation and scat¬ 
tering of neutralinos with SM particles helps us identify correlations between M p h and the SUSY 
spectrum. These correlations could be very helpful for constraining Tkd indirectly from current 
DM experiments. In particular, the region of m x o < 600 GeV could be potentially tested by the 
LHC, as commented in Appendix [Aj 


IV. IMPLICATIONS FOR DIRECT DETECTION 

Direct detection experiments of DM look for energy deposition in the underground detector 
caused by scattering interactions between target nuclei and WIMPs around us. The measurement 
or the bound on this cross section has direct consequences on the value of the T^, assuming that 
the processes involved in the last scattering are the same as the ones mediated the scattering of 
the DM with the detectors. 


Ni 2 W a + N 13 H 1 + 7 V 14 R 2 ■ The gaugino fraction is defined by Z g = |Yn| 2 + |A^i 2 1 2 (see [1001 for details). 
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Reference m analyzed correlations between the mass of protohalos, M p h (as well as the tem¬ 
perature of kinetic decoupling, T^d) and the SD and SI scattering cross sections. Such a correlation 
appears when the mass of squarks is assumed to be large ( rriq ~ 5-10 TeV)p*| and the dominant 
process for the scattering is mediated by a Z boson. 



o 


O 



_48 -46 -44 -42 


-52 -50 -48 -46 -44 


log w (a^ D (cm 2 )) 


logi^Gp 1 (cm 2 )) 


FIG. 3. Most probable regions on the scattering cross section-protohalo mass plane. The left and right 
panels show the correlation with the spin-dependent and spin-independent cross sections, a^ D and Up 7 , 
respectively. 


In Fig. [3j we show the most probable region on the plane of the protohalo mass, M p h, and the 
SD and the SI cross sections computed at tree level. Contrary to Fig. [l} the probability regions do 
not have disconnected parts, but they include both Higgsino-like (at ~1 TeV) and Wino-like (at 
~3 TeV) neutralinos. In both cases the dominant scattering process is mediated by the Z-boson. 
We see how the expected improvement on the SI sensitivity by, e.g., XenonlT [102j and LUX- 
Zeppelin experiment (LZ) [103| , will reduce the most probable range of the minimal subhalo mass 
down to below ~10~ 9 Mq, while the expected SD sensitivity provides weaker constraints. 

Since in the analysis we have included the XENON 100 limits as constraints on the WIMP- 
nucleon scattering cross section, we see in the right panel of Fig. [ 3 ] that the region around a ~ 
2 X 10 -44 cm 2 is strongly penalized. The current LUX bound is more stringent on the spin- 
independent sensitivity, giving an upper bound of a^ 1 ~ 10~ 44 for a 1 TeV neutralino [1104) . 
although we did not include it in our analysis. Including the LUX bound, therefore, would affect 
the very right part of the right panel of Fig. [ 3 ] (and also Fig. [ 5 ] shown below). However, since the 
regions affected are tiny, it would not affect our conclusions. 

Figure [4] shows points that reproduce the experimental constraint at 2 a confidence level for the 
minimal protohalo mass versus the tree level SD cross section plane. The right panel shows the 
case where the lightest first or second generation of sfermions is at least nine times heavier that the 
lightest neutralino, A(m^- —m^o) > 0.8. The thin yellow line corresponds to ~1 TeV Higgsino-like 
neutralino, while the thin red line to ~3 TeV Wino-like neutralino. In these two cases the Z-boson 
mediates both scattering processes. The rest of the points correspond to the Bino-like neutralino 
where, instead of a line, we get scattered points with 100 GeV < m v 0 < 1 TeV. We remind that for 

A-l 

16 The authors of Ref. Em used the squark mass to show the effect of light sfermions in the correlation, but clarify 
that when the correlation is broken, the relative contribution from squark, especially the slepton exchange in the 
kinetic decoupling process, increases. 
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A (mjq — m x o) < 0.8 


A (mj — m x <>) < 0.8, A (rriq — m ^o) > 0.8 


A (mjq - w^o) > 0.8 
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FIG. 4. Points that reproduce all the experimental observables at 2er confidence level in the SD cross section 
<jp D versus protohalo mass M p h plane. The neutralino mass is indicated with colors, as shown in the color 
bar. The three panels separate the points in three groups: light squarks and sleptons (left panel), light 
sleptons and decoupled squarks (central panel), and decoupled squarks and sleptons (right panel). 


the Bino-like case the annihilation cross section and, therefore, the relic density can be adjusted 
varying the neutralino mass and its mass splitting with the lightest (Wino-like) chargino. On 
the other hand, even if it is ten times heavier than the lightest neutralino, sleptons mediate the 
dominant scattering processes that set T^d for most of the points. The size of the contribution 
of processes, mediated by the Z-boson, depends on how large the Higgsino component of the 
neutralino is. However, the Higgsino component of a Bino-like neutralino is highly constrained by 
SI cross sections bounds. Nevertheless, as we comment in subsection III C[ there are some blind 
spots for SI cross sections. For those points the Z-boson gives an important contribution to the 
scattering cross section. 

Regarding for the Bino-like region, the dominant process is mediated by the Z-boson p^j 
Besides the dominant scattering processes for Tkd and er® D are different, there is an apparent 
correlation between the two quantities for a fixed neutralino mass. We have checked the behavior 
of the correlation for a specific values of m x o, finding that spreads around one order for a given 


value of o|i D . 

Another characteristic of the Bino-like case, assuming sfermions are heavy, is that the minimal 
protohalo mass is the one allowed by the free-streaming. On the contrary, for the Higgsino-like 
and Wino-like cases, the minimal protohalo mass is the one allowed by the damping scale set by 
acoustic oscillation. 

Central panel of Fig. [4] shows the case where the lightest slepton has a mass smaller than ~10 
times the lightest neutralino mass. As expected, the Wino-like and Bino-like regions spread to 
larger protohalo masses^ Left panel of Fig. [I] shows the case where the lightest sleptons and 


17 Squarks are typically heavier than sleptons when parameterizing the model at gauge coupling unification scale. 
Therefore, imposing the condition A(m;- — m x °) > 0.8 implies that squarks are typically much heavier than ten 
times the mass of the lightest neutralino. 

18 Winos and Binos have strong SD interaction since diagrams where the incoming and outgoing fermions have the 
same helicity are allowed. On the other hand, the diagrams where the incoming and outgoing fermions have opposite 
helicities are spin-independent, requiring a qqH vertex to yield the helicity flip, which is Yukawa suppressed. For 
a review, see Ref (105] . 















16 


squarks are smaller than ~10 times the lightest neutralino. Here, squarks are light enough to give 
important contributions to the scattering with the nucleus, spreading the points to larger values 
oia s n D . 


A (mjq — m < 0.8 


A (mj — m^)) < 0.8, A (m$ - m x o) > 0.8 


A (mjj - m x o) > 0.8 
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FIG. 5. Same as Fig. [l] for the the SI cross section, a SI . 

Figure [5] shows points in the minimal protohalo mass versus tree level SI cross section plane. 
The main contribution to the SI cross section comes from the Higgs exchange, requiring a non- 
negligible Higgino and Wino/Bino coupling (since Higgs couplings through neutralinos are HHB 
and HHW). On the other hand, the total neutralino-SM scattering, and therefore Tkd and M p h, 
are dominated by SD interactions. A a consequence, right panel of Fig. [5] shows the correlation 
between the ZxiXi and Hx\Xii f° r the Higgsino and Wino case. Central and left panels show the 
effect of sleptons and squarks in the scattering processes. 

Figures |4] and [5] show the expected sensitivity by XenonlT and LZ assuming the neutralino 
mass is ~l TeV. For a neutralino of ~T00 GeV, the expected sensitivity is around one order of 
magnitude stronger. 

As we commented in the previous section, assuming universality and unification of the squark 
masses and slepton masses, we impose a particular mass hierarchy: t\ is the lightest squark and f\ 
is the lightest slepton. Without this assumption, first and second generation of sfermions can be 
lighter and change the phenomenology for direct detection experiments and colliders. In the case 
that the lightest neutralino is gaugino-like and the first and second generation of squarks are the 
lighter sfermions, Tkd will be completely correlated with the neutralino-nucleon scattering cross 
section. Still, if they are not the lightest ones but they are significantly lighter than in our analysis 
the neutralino-nucleon scattering cross section could increase, getting values close to the actual 
limits. In addition, if they are lighter or close in mass to the first and second generation sleptons, 
the scattering of the neutralino with SM particles in the early universe could also increase. Those 
points will most likely populate the top-right corner of the left panel of Figs. [4] and [5j 

Notice that, as we mentioned above, for the computation of SI cross sections we have adopted 
fTs = 0.36 for the contribution of the strange quark to the nucleon form factors m, derived 
experimentally from measurements of the pion-nucleon sigma term. However this value is con¬ 
siderably larger than determinations obtained from lattice QCD, fx 3 = 0.043 ± 0.011 [ 106] . The 
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discrepacy between the two values and its impact in the SI cross section is studied in more detail 
in Refs. mm- 

Finally, we remind that the scattering cross-sections considered in this work were computed at 
tree level. In the cases where neutralino approaches to a pure state (Bino, Wino or Higgsino), this 
approximation may not give a reliable result. In particular, in the case of the Wino-neutralino, 
one loop corrections give the dominant contributions (see, e.g. Ref. [Ml ED]). 


V. IMPLICATIONS FOR INDIRECT DETECTION 

One of the most reliable methods to model the non-linear evolution of DM is numerical simu¬ 
lation, although it is limited by mass resolution. In fact, the minimum self-bound mass (M nun ) of 
DM halos is expected to be many orders of magnitude below the resolution of current simulations. 
Through numerical simulations such as Acquarius EH, we can obtain information on the subhalo 
hierarchy, although its resolution mass limit 10 4 or 10 5 Mq is far from the predicted protohalo 
mass shown in Sec. EO 

Here we investigate the impact of different values of M m j n on the gamma-ray luminosity due to 
DM annihilation, and compute a boost factor of a given halo of mass M due to the substructure 
inside it, by integrating the subhalo annihilation luminosities from the protohalo mass we have 
found, M p h, up to the mass of sizable fraction of the host halo M max . The total luminosity of the 
DM halo due to annihilation is proportional to: 

f* ATiax .//) 

/ Mpl dM m L - ‘M- < 7) 

where dn/dM is the subhalo mass function, i.e. the subhalo number density per unit mass range. 
Numerical simulations find that the differential subhalo mass function follows a power law dn/dM oc 
M - - 5 , with /3 ~ 1.9 or [3 ~ 2 (see, e.g. imiEa). We adopt a M 2 subhalo mass spectrum as our 
fiducial subhalo model. 

We assume that each individual DM subhalo is described by a Navarro-Frenk-White (NFW) 
density profile [ 114j : 


Psh {r/r s )(l + r/r s ) 2 ’ ^ 

where p s and r s are the characteristic density and radius, respectively. L s h(M) is defined as the 
luminosity of each subhalo in the host halo, which depends on the volume integral of the subhalo 
density squared, and is given by: 


L sh (M) = J dV sh p 2 sh oc p 2 s r 3 s . 


(9) 


Following the same approach of Ref. EH, we parameterize the scaling relation between the 
gamma-ray luminosity and subhalo mass as: 


L sh (M) oc L 0 x 


( M 

l io 4 m q 

l M 

V 10 4 Mq 


0.77 

7 


M > 10 4 Mq 
M < 10 4 Mq , 


( 10 ) 


where above the simulation resolution of ~10 4 Mq, the luminosity versus subhalos mass scales 
as L oc M 0 77 , while below the resolution we assume 7 < 1. Here Lq encodes all the particle 
















18 


physics, i.e., Lq oc ( av)/m 2 0 , where (av) is the velocity-averaged annihilation cross section times 

^1 

the relative velocity| 19 | 

In order to obtain the scaling behavior of L s h oc M 0,7 , we adopted scaling relations among 
several quantities found in the Aquarius numerical simulation. Since each subhalo is described by 
a NFW density profile, we related the maximum rotation velocity of the subhalo, V ma , x , and the 
radius at which the rotation curve reaches this maximum, r max , with the characteristic density and 
radius, p s and r s , to obtain them as a function of the subhalo mass M. These empirical relations 
between (U ma x, r m ax) and ( p s , r s ), however, lose validity in mass regions below the resolution limit 
of the simulation. For this reason we split Eq. (10) in two terms, above and below the resolution 
(10 4 Mq), where in the latter we put 7 as a phenomenological parameter describing the scaling 
behavior. 


The luminosity in Eq. 0 can be then written as: 
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Assuming that the first term dominates, the luminosity is, thus, a function of the protohalo mass: 


L(M ph ) 
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For comparison, we define a reference value for such a luminosity, L re f, as: 

M ref 
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m 0 

x? 
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( 12 ) 
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For values of these reference parameters, we adopt (av) re f = 3 x 10 -26 cm 3 s , M re f = 10~ 6 Mq, 
and 7 = 0 . 8 . 

The left panel of Fig. [ 6 ] shows the two-dimensional joint posterior PDF for the protohalo mass 
M p h and (av) with 68 % and 95% credible contours. These most probable regions fall in a mass 
range between 10 - ' and 10 ~ 12 Mq, and (av) = 10 _ 26 -10 -24 cm 3 s _1 . The region with higher 
probability density again corresponds to a Higgsino DM candidate with the annihilation cross 
section close to the canonical value 10 _ 26 cm 3 s _1 , while the second region corresponds to a DM 
Wino candidate with much larger annihilation cross section ~1CU 24 cm 3 s . In the right panel we 
show the ratio of the luminosity over the reference one L = L/L re f , versus the DM mass, m x 0 . 
We also analyzed the change in the boost by varying the 7 -paranreter in a range between 0.5 and 
0.9, we only show the case 7 = 0.8, and found that L always got largely boosted by decreasing 7 . 
This behavior depends on the normalization made on the protohalo mass, M p h, since it has been 
normalized to the limit of the numerical simulation (10 4 Mq). 

Figure [ 7 ] shows the boost factor, L = L/L re f, for points that reproduce all the experimental 
observables within 2 a confidence level. Right panel shows points which refer to a Higgisino-like and 
Wino-like neutralinos, while the left panel shows points where the neutralino is mostly Bino-like. 
Bino-like neutralinos have very small (av) in the limit of zero velocity. Co-annihilations, which play 
a very important role in the efficient annihilation in the early Universe, are not present anymore; 
this is the reason for which we have a very small boost of the luminosity. 

Finally we comment that although not included in this work, Fermi and HESS bounds on 
m. x o-(av) plane strongly constrain the wino-like region, excluding the region around 2.4 TeV; see 
Refs. pTmiTT9] . 
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I 11 the considered MSSM, for almost all the data points, we find that the annihilation cross section, (av), is almost 
independent of velocity, (av) & (av)o- 
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logio[(av)(cm 3 s 1 )] logio(myo(GeV)) 


FIG. 6 . The two dimensonal joint posterior probability density function for the protohalo mass, M p h, versus 
the velocity-averaged annihilation cross section times the relative velocity, (crv) (left panel), and for L, 
obtained by using 7 = 0.8, versus the Dark Matter particle mass, m x 0 (right panel). For both panels, the 
region with higher probability density corresponds to a Higgsino DM candidate; in the second region the 
DM candidate is a Wino. Left panel shows that, in the Higgsino case, the protohalo mass M p h is lower 
than the reference one, while (av) does not deviate from (av) le { ~ 10 -26 cm 3 s _1 . Right panel shows that, 
in the Wino case, the protohalo mass M p h is even lighter and (av) is two orders of magnitude larger than 
(av) le { ~ 10 ~ 26 cm 3 s _1 ; thus, there is a substantial enhancement of L. 



FIG. 7. The mass of the lightest neutralino versus the boost factor, L = L/L re f, for points that reproduce 
all the experimental observables within 2er confidence level. Left panel shows points which refer to a Bino¬ 
fraction (Nn) larger than 0.8. Right panel shows point with a Bino-fraction smaller than 0.8 



VI. SUMMARY AND CONCLUSIONS 

In this work we have studied how the kinetic decoupling of dark matter could improve our 
knowledge of the properties of the dark matter protohalos within a supersymmetric model, i.e. the 
Minimal Supersymmetric extension of the Standard Model. Such a model is the well-motivated 
extension of the Standard Model at the electroweak scale. At first, it was introduced to solve 
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the hierarchy problem of the Standard Model, but it revealed to have many other interesting 
characteristics. In particular, it contains a tempting particle dark matter candidate, the lightest 
neutralino. If such a neutralino is the lightest supersymmetric particle and the quantum number 
R-parity is preserved, it is stable, yielding to a thermal abundance as that indicated by the observed 
dark matter density. 

In our analysis we do a forecast on the mass of the protohalos within a supersymmetric frame¬ 
work realized with 9 independent parameters. We performed any analysis in the light of the latest 
data coming from particle physics experiments, as well as the relic density constraints. Among 
them, the most important observables involved in the analysis, which give a relevant impact on 
our results, are the mass of the Higgs and the relic density. 

1. The kinetic-decoupling temperature and, thus, the minimal protohalo mass result to be not 
well constrained for WIMPs, since the interactions involved in the annihilation of neutralinos, 
that are constrained by the relic density, are not necessarily those which participate in the 
scattering of neutralinos with first and second generation of fermions. In a supersymmetric 
framework, the minimal protohalo mass is typically 1CU 6 Mq, assuming a Bino-neutralino 
annihilating through sfermions with a mass of around twice the neutralino mass. This 
resulted in a possible option to get a well tempered neutralino. In addition, this possibility 
has been well motivated by constrained scenarios like CMSSM, affirming that when the 
neutralino is mostly Bino, it efficiently annihilates through sfermions in the early Universe, 
giving the correct relic density. Nevertheless, it was in tension with the experimental data 
within the CMSSM, especially after the first run of the LHC, where a considerable part of 
this region was excluded. 

2. Using a Bayesian framework, we showed that the most probable neutralino mass regions 
satisfying both the Higgs mass and the relic density contraints, are those with the lightest 
supersymmetric neutralino mass around 1 TeV and 3 TeV, that correspond to Higgsino- 
like and Wino-like neutralino, respectively. We mentioned that, concerning the Higgsino-like 
neutralino, the annihilation cross section is driven by its Higgsino component, while for Wino- 
like neutralino, the annihilation cross section is mainly driven by its Wino component. We 
also discussed that the part of the region closer to ~ 2.4 TeV gets important contributions 
from the neutralino-stau co-annihilation, reducing both the effective annihilation rate of 
neutralinos in the early Universe and the value of the neutralino mass, in order to obtain 
the correct relic density. 

3. We commented that in the case of Wino-like or Higgsino-like neutralinos the annihilation 
products are gauge bosons, whose interactions involve different couplings with respect to the 
ones of the neutralino-fermion scattering. For that reason kinetic decoupling temperature, 
Tkd, exhibits a considerable range of variation, that reflects, in turn, to a protohalo mass 
range of M p h ~ 10 _12 -10 -7 M 0 . 

4. We also discussed the Bino like neutralino with masses smaller than ~1 TeV, where a quasi¬ 
degenerated sfermion or chargino, or a light sfermion are necessary to get the correct dark 
matter abundance. Sleptons give the most important contribution for the kinetic decoupling 
temperature and therefore to the protohalo mass, setting the range M p h 10 -11 -10 -4 M q . 

5. Kinetic decoupling of dark matter, involving elastic scattering of a dark matter particle 
with Standard Model particles in the early Universe, reveals a relevant process for dark 
matter direct detection searches. In our analysis, we showed that the regions where the 
probability is higher the correlation between the protohalo mass and experimental signa¬ 
tures permits to put constraints on the protohalo mass. We depicted how improvements 
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on the spin-independent sensitivity might reduce the most probable range of the protohalo 
mass between ~1CR 9 Mq and ~10 - ' Mq, while constraints associated to the expected spin- 
dependent sensitivity are weaker. To give this conclusion we computed scattering cross 
sections at tree-level. However, specially in the Wino-like neutralino case, loop corrections 
should be considered since the tree level coupling vanishes when approaching the pure Wino 
case. 

6 . We discussed, as well, how the interplay among both spin-dependent and spin-independent 
scattering procesess, strongly depends on the neutralino composition. For both Higgsino- 
like and Wino-like cases, the spin-dependent scattering between Higgsino and fermions is 
mediated dominantly by the Z boson at tree level, while for the spin-independent scattering, 
the interaction is mediated by the Higgs boson. Regarding the Higgsino neutralino, we 
commented that the spin-independent interaction gives a nonzero tree-level contribution as 
long as gauginos are not decoupled, a non-negligible Bino or Wino component is necessary 
to have a non-negligible coupling with the Higgs. On the other hand, for the Wino-like 
neutralino the requirement of a non-negligible component of Higgsino is indispensable to 
have a tree-level contribution to both scattering processes if sfermions are decoupled. 

7. Depending on the nature of neutralino, the value of the annihilation cross section, (av) v ^o, 
changes by different orders of magnitude. We presented that the annihilation cross sec¬ 
tion, (av), in the Higgsino case does not deviate from the canonical cross section, (av) ~ 
10“ 26 cm 3 s -1 . On the other hand, in the Wino case non-pertubative effect is important, 
(av) increases up to two orders of magnitude. And it is much smaller in the Bino-like case, 
where co-annihilations with sfermions played a crucial role to fix the correct abundance. 

8 . Another way to look for dark matter is through indirect detection methods, which consist to 
detect, indirectly, the lightest supersymmetric particle through annihilation processes where 
Standard Model particles, including gamma-ray photons, are produced. Since the luminosity 
of each subhalo in the host halo due to the dark matter annihilation processes depends on 
the volume integral of the subhalo density squared, smaller and denser substructures provide 
an enhancement of the luminosity. In this work, we showed for both neutralino Higgsino- 
like and Wino-like cases how the boost of the luminosity due to dark matter annihilation 
increases, depending on the protohalo mass. We discussed that in the Higgsino case, there 
is no a significant enhancement of the luminosity: the protohalo mass is lower than the 
standard value often used in the literature of ~10 -6 Mq, while (av) does not deviate from 
(av) ~ 10 -26 cm 3 s _1 . In the Wino case, a substantial enhancement of the luminosity is 
seen: the protohalo mass reaches lighter values, and (av) is two orders of magnitude larger. 
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Appendix A: Implications for collider searches 


As commented in Sec. 


me 


points with m v o smaller than ~1 TeV are Bino-like and require a 

i—i . ^1 

light enough next-to-lightest sparticle, in order to guarantee an efficient annihilation in the early 
Universe. Based on the characteristics of the next-to-lightest sparticle, we are going to comment 
the potential LHC signatures. 

For neutralinos lighter that 500 GeV there are two regions, in addition to Z/h/A resonances. The 
first one has X\ close in mass to Xi- A light Wino-like chargino which annihilates and co-annihilates 
in the early Universe is required , and is represented by points with 5 GeV < m x o — m x o < 40 GeV in 
the left panel of Fig. [8j In this region Xi is dominantly Bino and xf and X 2 are dominantly Winos. 
The Bino and Wino mass, M\ and M 2 , are close to the values where the tree level decay of X 2 t° 
zw Xi is suppressed, and the branching ratio to jXi squires a large value, as discussed in detail in 
Refs. il20 . H2~l : . Right panel of Fig. [8] shows that some of the points can have a dominant X 2 7Xi 
decay, giving a characteristic signature at collider. Moreover, the decay channel II —> lx® ^7Xi 
becomes relevant. Although the photon produced in the X 2 an d II decays is very soft, it could give 
a clear signature at collider in the boosted regime. Keep in mind that a potential measurement of 
sleptons will directly constraint the prediction for the protohalo mass for Bino-like neutralino. 



> 

D 
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FIG. 8. Left panel shows the Xi mass and X 2 mass plane. The colors show the mass of the lightest slepton 
{Sl.r, ^l,r)- Right panel shows the branching ratio of X'i to photons as a function of X 2 mass for points 
with mi L > m x 0 . Red, blue and green points correspond to m 7 y > 1 TeV, 500 GeV < m ; y <1 TeV and 
rrij < 500 GeV respectively. 


The second region corresponds to stau co-annihilation, where f\ and \'5 are very close in mass. 
In the left panel of Fig. |8j the points outside 5 GeV < m x 0 — m x 0 < 40 GeV correspond to this 
region. Notice that, as a consequence of universality conditions of slepton soft masses, the first and 
second generation of sleptons is relatively close in mass to the lightest stau and, therefore, to the 
lightest neutralino. The authors of Ref. 1122 discuss the status of this region after the first run of 
the LHC in the framework of the constrained MSSM (CMSSM), and project the likely sensitivity 
of the LHC searches in Run 2 at 14 TeV center of mass energy and 300/fb of integrated luminosity, 
concluding that the entirely CMSSM co-annihilation strip will be tested. 

For m x 0 > 500 GeV new regions arise. Stop co-annihilations, and neutralino annihilations are 
mediated by sfermions. In this neutralino mass range the production of colored particles is the 
most promising. In Refs. {123H125] . is studied the stop co-annihilation region, not only by direct 
stop production but also by gluino production, where direct stop productions constraint fight stops 
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{ rn t 1 1$ 400 GeV); for heavier stops gluino, the production seems to be more promising. On the 
other hand, the region where neutralino annihilation is mediated by squarks is directly constrained 
by limits on squarks masses. 
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